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Abstract 

We establish a new formula for the heat kernel on regular trees in terms of classical J-Bessel 
functions. Although the formula is explicit, and a proof is given through direct computation, 
we also provide a conceptual viewpoint using the horocyclic transform on regular trees. From 
periodization, we then obtain a heat kernel expression on any regular graph. From spectral 
theory, one has another expression for the heat kernel as an integral transform of the spectral 
measure. By equating these two formulas and taking a certain integral transform, we obtain 
several generalized versions of the determinant formula for the Ihara zeta function associated 
to finite or infinite regular graphs. Our approach to the Ihara zeta function and determinant 
formula through heat kernel analysis follows a similar methodology which exists for quotients 
of rank one symmetric spaces. 

1 Introduction 

Let X be a (<7+l)-regular graph for an integer q > 0. There is an associated heat kernel Kx {t, xq, x) 
corresponding to the Laplacian formed by considering the adjacency matrix on X . Wc show that 
in a natural way the building blocks of Kx are the functions 

q -r/2 e -(^)t Ir(2Vq - t ), 

where r £ Z>o , time t £ R>o, and I r is the classical /-Bessel function of order r. The expression 
for the heat kernel on X comes from a new formula for the heat kernel on regular trees (Proposition 
13. 1[) which we prove in this article. Our expression is quite different from a previous formula due 
to F. Chung and S.-T. Yau jCY99j . which we describe in subsection 13. 1[ Equations (|5al I5b[) . If we 
write the functions in the above stated building block as 



<1 



-r/2 . e -(^-i) 2 « . e -^t Ir {2^q t ), 



then there is moreover a near-perfect analogy with the building blocks of typical heat kernel 
expressions on Ricmannian symmetric spaces, which have the form 

F(r) ■ er at 




for certain constants a, d and function F that allow further interpretations; we refer to the survey 
article jJLOlj and the references therein for further discussion. 

More precisely, we prove the following result. 
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by all authors. We thank Pierre de la Harpe for a number of corrections. 
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Theorem 1.1. The heat kernel on a (q + l)-regular graph X is given by 



K x (t,x ,x) = e-^ +1 )' Y, b m (x)q- m/2 I m (2Vqt), 

m=0 

where I m is the I-Bessel function of order m, b m (x) = c m (x) — (q — l)(c m _2(2;) + Cm-i{x) + •■■) 
and c m (x) is the number of geodesies from a fixed base point xo to x of length m > 0. 

To be specific, we define a geodesic in a graph to be a path without back-tracking. The terminology 
is consistent with concepts from Riemannian geometry where a geodesic is a path which locally 
is distance minimizing. Moreover, a closed geodesic is a closed path without back-tracking or 
tails. We defer until section [2] for more details on these definitions and for a precise axiomatic 
characterization of the heat kernel. In view of a combinatorial observation, Proposition 12.11 we 
may formulate the following result. 

Corollary 1.2. In addition to the stated assumptions, suppose that X is vertex transitive. Let 
denote the number of closed geodesies of length m in X with base point xq . Then 

oo 

K x (t,x ,x ) = K q+1 (t,x ,x ) + e-^ t N^q~ m / 2 I m (2^t) 

m— 1 

where K q+ \ denotes the heat kernel of the (q + \)-regular tree. 

For finite, not necessarily vertex transitive, graphs X, similar formulas were previously proved by 
Ahumada }Ahu87| . Terras- Wallace |TW03j and Mnev |Mne07) . Note that our formula holds for 
infinite graphs as well, and is therefore more general than the finite graph case. With our methods, 
one can also deduce the formula for the finite non- vertex transitive case using Proposition 12.21 

There is a second expression for the heat kernel coming from spectral considerations. Equating the 
two expressions for the heat kernel, as in known approaches to the Poisson summation formula or 
the Selberg trace formula, one obtains an identity which is a type of theta inversion formula. From 
the identity, we will apply a certain integral transform, which amounts to a Laplace transform with 
a change of variables from which we obtain the logarithmic derivative of the Ihara zeta function. 
This procedure is motivated by McKean |McK72| in his approach to the Selberg zeta function and 
was axiomatized in |JL01j to abstract settings. In the end, one obtains determinantal formulas 
for Ihara zeta-likc functions. In a special case, we recover the standard formula stemming from 
Ihara's work |Iha66j . which in turn is generalized in jHas891 IBas92[ IFZ991 IKS001 ITS07] for finite 
graphs; see subsection 14.31 

We now describe one sample outcome which comes from the above described sequence of calcu- 
lations. Let X be a vertex transitive (q + l)-regular graph. We define the associated Ihara zeta 
function of X by 

( °° M° 1 
CxW=ex P £-f« m , 
Lm=l J 

where is the number of closed geodesies of length m starting at a fixed vertex xq. For finite 
graphs, the classical Ihara zeta function is just our Ihara zeta function raised to the power equaling 
the number of vertices. 

Theorem 1.3. Let X be a vertex transitive (q + \)-regular graph with spectral measure [i for the 
Laplacian. Then 

(xiu)- 1 = (l-u 2 ) (9 - 1)/2 cxp (J \og(l-{q+l-\)u + qu 2 )dii(\)^ . 

Again, we defer to section [3] for the definitions of the Laplacian and spectral measure. There are 
some papers in the literature defining Ihara zeta functio ns for infinite graphs, in particular Clair 
and Mokhtari-Sharghi |CMS01| . Grigorchuk-Zuk |GZ04| and Guido, Isola and Lapidus |GIL08j . 
Their definitions are at least a priori somewhat different in that they typically look at approxi- 
mations by finite graphs and also using von Neumann trace for group operator algebras. In many 
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cases it coincides with our definition and thereby the formula in Theorem If .31 can be recovered in 
those references. 

A numbe r of interesting examples of Ihara zeta functions for infinite Cayley graphs can be found 



in [GZ04j. Additionally, we refer to the articles [Sun86j . |Sun94] and |Sun08| in which the author 
provides a fascinating discussion in which the spectral and zeta function analysis on graphs is 
compared to similar studies in spectral theory on symmetric spaces and zeta functions from number 
theory. 

In summary, we have given a new expression for the heat kernel associated to any regular graph. 
One can quickly deduce the Ihara determinant formula and a number of interesting extensions, 
not the least of which is to infinite transitive graphs. As is well-known, one main application of 
such formulas is to the study of counting closed geodesies. For us, it is also significant that our 
analysis provides yet another instance of when the heat kernel yields zeta functions together with 
their main functional relation, just as in the case of Riemann, Selberg, and beyond. In particular, 
the present paper can be viewed in the context of the last section of [JLOlj . 



2 Preliminaries 

2.1 Graphs 

We follow the definitions in Serre's book |Ser03j . A graph X consists of a set VX which are called 
vertices, a set EX which are called edges, and two maps 

EX^VXxVX, y h+ (o(y),t(y)) 

and 

EX EX, y^y 
such that for each y £ EX we have that y = y, y ^ y and o(y) = t(y). 

The vertices o(y) and t(y) are the extremities of the edge y. Two vertices are adjacent if they are 
extremities of an edge. The degree of a vertex x is 

degx = Card{y G EX : o(y) = x}. 

A graph is d-regular if each vertex has degree d. 

There is an obvious notion of morphism. Let PATH n denote the graph with vertices 0, 1, 2, ...,n 
and (half of) the edges arc given by [i, i + 1], i = 0, n — 1. A path (of length n) is a morphism c 
from PATH n into the graph. The sequences of edges yi = c([i, i + I]) (such that t{yi) = o{yi + \)) 
determines the path. In particular a path is oriented. There is a backtracking if for some i that 
Vi+i = yl and there is a tail if yo = y n _i- A path is closed if c(0) = c(n). A geodesic is a path 
without backtracking. A geodesic loop (or circuit in Serre's terminology) is a closed path that 
is a geodesic. A closed geodesic is a closed path with no tail and without backtracking. (This 
is analogy with Riemannian geometry where a closed geodesic, as opposed to a geodesic loop, is 
required to be smooth also at the start/end point). The path of length zero counts as a closed 
geodesic and, therefore, is a geodesic loop. Additionally, every closed path with one edge counts 
as a closed geodesic. Any length two geodesic loop is also a closed geodesic, but the closed path 
y.y is neither. 

A prime geodesic is an equivalence class of closed geodesies [c], where the equivalence class is 
forgetting the starting point and which is primitive in the sense that it is not a power of another 
closed geodesic. The latter means by definition that there is no closed geodesic d and integer n > 1 
such that [c] = [d n ], which says in words that c is not just a geodesic that traverses another one 
n number of times. 

An orientation is a subset EX + of edges such that EX is the disjoint union of EX + and EX + . 
With the data of a graph one can associate a geometric realization: start with the discrete topology, 
take VX x [0, 1] and make identification based on the maps o and t. 

A tree is a connected nonempty graph without geodesic loops. 
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We will in particular consider vertex transitive graphs, that means that there is a group of auto- 
morphisms which is transitive on the vertices. In particular such a graph is of course regular. A 
rich source of such graphs is provided by Cayley graphs of groups: Let G be a group and let S be 
a subset of G. We denote by X(G, S) the oriented graph having G as vertices and EX+ = G x S 
with o{g, s) = g and t(g, s) = gs for each edge (<?, s). 

Let X be a graph on which a group G acts. An inversion is a pair consisting of an clement g and 
an edge y such that gy — y. If G acts without inversions (which is the same as saying that there 
is an orientation of X preserved by G) we can define the quotient graph G\X in an obvious way; 
the respective edge and vertex sets are the corresponding quotients. (To get rid of inversions one 
may pass to a barycentric division.) As in topology we say that X is a regular covering of Y if 
there is a group which acts on X freely and without inversion with quotient Y. 
Let X be a (g + l)-regular graph. Then its universal covering is the (q + l)-regular tree, and the 
covering group acts freely on the tree without inversion. The covering group is a free group. 



2.2 Path counting in graphs 

We fix a base vertex xq in a graph X and define the following counting functions which will be 
used in this paper: 

• afe(x) is the number of paths of length k from xo to x, 

• Cfe(x) is the number of geodesies of length k from xq to x, 

• c° = Ck(xo) is the number of geodesic loops of length k starting at Xq, 

• Cfc the number of geodesic loops of length k, from some starting point with a distinct direction, 

• iV° is the number of closed geodesies of length k starting at xq, 

• Nk is the number of closed geodesies of length k, from some starting point with a distinct 
direction, 

• 7Tfc is the number of prime geodesies of length k. 

The sequences {ck}, {Nk} and {irk} only make sense for finite graphs, since if the graph is infinite, 
these values are typically infinite. For finite vertex transitive graphs, the sequences {N®} and 
{Nk} are related by the number of vertices, i.e. starting points. Specifically, if the graph X has n 
ver tices, t hen N® • n = Nk for all k. Also, Nk and TTk have a precise relationship, see e.g. [Terllj 
or 



GZ04j. Finally, we recall the conventions that eg = N$ = 1, c? = A^ and c\ = N$. 



Proposition 2.1. Let X be a transitive (q + l)-regular graph. Then for k > 3. following relation 
holds true 

N° = 4-(q-l)(4_ 2 +4_ i + ...) 
the last term being 4 or c 2 depending on the parity of k. 

Proof. This is similar to an argument in [Ser97\ A geodesic loop of length k > 3 which is not 
a closed geodesic has the form yi.z.y 1 where z is a geodesic loop of length k — 2. There are two 
possibilities, either z is a closed geodesic or not. If we fix z, then the number of possibilities for y\ 
is q — 1 in the first case and q in the second case. Since X is vertex transitive, we may freely change 
the starting point of any loop, namely z. With this in mind, we obtain the recursive relation that 

4 AT° = (q - l)N° k _2 + <l(4-2 N° k _ 2 ) = (4_ 2 N°_ 2 ) + (q l)c°_ 2 , 

which we can write as 

Nt-Nl_ 2 = cl-qcl_ 2 . 
Using that c\ = N® and c 2 = N 2 , the proposition follows by induction on k. □ 

With a proof similar to the one given in the above proposition, we obtain the following result. 
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Proposition 2.2. Let X be a finite (q + l)-regular graph. Then for k > 3, the following relation 
holds true 

N k = c fc - (q - l)(c fc -2 + c fc _ 4 + ...) 
the last term being c\ or C2 depending on the parity of k. 

2.3 The combinatorial Laplacian and heat kernel 

Given a (q + l)-regular graph X and function / on the vertices X, the Laplacian of /, written as 
A/, is the function of the vertices of X which is defined by the formula 

Af(x) = (q+l)f(x) - £ f(t(e)). 

e s.t. o(e)—x 

The Laplacian is a semi-positive, bounded self-adjoint operator on L 2 (VX). For a finite graph 
with N vertices we label the eigenvalues of A as follows: = Ao < Ai < ... < Ajv_i < 2(q + 1). 
The heat kernel Kx(t, x, y) : R>o x!xI->IonIis the solution of 

d 

AK x (t,xo,x) + —Kx{t,x o ,x)=0 (la) 
ot 

K x (0,x ,x) = { 1 [ix = Xo (lb) 
I otherwise. 

Sometimes we write Kx(t,xo,x) = Kx(t,x) when a base point x$ is understood. 
Whenever X is a countable graph with bounded vertex degree, the heat kernel of X exists and 
is unique among bounded functions QDM061 lDod06) ). Let X be a regular covering of Y via the 
map 7r and group L, then it is formally immediate that 

K Y (t,yo,y)= ^2 K x (t,x ,x) 

where ya = ir(xo); alternatively, with a chosen x such that y = tt(x), we have that 

Ky{t,y ,y) = ^Kxit^o^x). 

One can show that the heat kernel decays sufficiently rapidly so that the above equalities are not 
only formal, but indeed are convergent series. We refer to |CJK10| IDM06] as well as the heat 
kernel formula for trees in Proposition 13.11 and the bounds in subsection 12.51 

Consider the numbers a n (x) defined by 

e (i+^K x (t,x) = y / a n (xy- 1 

then it is well-known and simple to see that a n (x) is the number of paths from Xq to x as defined 
above. 

2.4 The i-Bessel function 

Classically, the I-Bessel function I x (t) is defined as a certain solution to the differential equation 

9 d 2 w dw . 9 9 , 

For integer values of x, it is immediately shown that I x = I^ x and, for positive integer values of 
x, we have the series representation 

Kt)-T {t/2)2n+x (2) 

71 — ' 
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as well as the integral representation 



I x (t) = -J e tcos(e) cos(8x)d0. (3) 
o 

The mathematical literature contains a vast number of articles and monographs which study the 
many fascinating properties and manifestations of the /-Bessel functions, as well as other Bessel 
functions. The connection with the discrete heat equation comes from the basic relation 

I x+1 (t)+I x - 1 (t) = 2^I x (t), (4) 

which easily can be derived from the integral representation and trigonometric identities. This 
relation will be used in the proof of Proposition 13.11 

2.5 Universal bounds for the /-Bessel function 

We have the following uniform bounds from ICJKlOj which used jPa!991 . For any t > and integer 
x > 0, we have that 

^- ^ (rb)" 2 =( i+ Tr /2 - 

As stated, the above bound is enough to show that the periodization procedure in our setting 
gives rise to convergent sum expressions for the heat kernel. 

2.6 An integral transform of I-Bessel 

For integers n and s€C with Rc(,s) ^ we have e.g. from |OB73j . that 

/ e- st e- l I n {t)dt = ± '—. 

We will consider the transform, essentially the Laplace transform, 

poo 

Gf{u) = (u~ 2 -q) e - [qu+1/u)t e (q+1)t f{t)dt. 
Jo 

In view of the above formula, applying the transform to the heat kernel building block, we get 

G(e-(« +1 V fe/2 4(2Vgi)) (u)=u k - 1 

for k > and u > 0. 

3 Heat kernels on regular graphs 

3.1 A heat kernel expression for regular trees 

Let X be the (q + l)-regular tree and Xo G X a base point. From its characterizing properties (jlajl 
and (|lbp. it is immediate to show that the heat kernel on a graph is invariant with respect to any 
graph automorphism g: 

K(t,gx ,gx) = K(t,x ,x). 

In particular, we have that the heat kernel K (t, xq, x) on the tree X is radial, that is it depends only 
on r = d(xo,x). Therefore we can write the heat kernel as K(t,r). Expressions for K(t,r) were 
established by Bednarchak jBed97j . Ch ung and Yau |CY99j . Cowling, Mcda and Setti jCMSOOj . 
as well as Horton, Newland and Terras HNT06 . In the physics literature regular trees are called 
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Bethe lattices. As stated, one of the main results of this paper is a formula for the heat kernel on 
X, which we consider to be new since we did not find the expression in either the mathematical 
or physics literature. 

As an example of a known expression, Chung and Yau |CY99j prove that in the radial coordinate 
r, the heat kernel of the (q + l)-rcgular tree is given by 

. , 2e~( q+1 ' )t f* exp (2t^/qcosu) smu(qsin(r + l)u — sin(r — l)it) , , 

K(t,r) = TT-r / — 7 7T5 : 5 du (5a) 

irq r / 2 - 1 J (<7 + l) -4gcos 2 u 



for r > 0, and 



2q(q + l)e-(9 +1 )* r exp (2t^qcosu) sin 2 ■ 



K(t,o) = ' ~'" / T - "~ a ' T du - ( 5b ) 

v ' ir J {q + l) 2 - 4q cos 2 u y ' 

The formula which we prove is given in the following proposition. 

Proposition 3.1. The heat kernel of the {q + 1) -regular tree is given in the radial coordinate r > 
as 

oo 

K(t,r) =q- r / 2 e-^ t I r (2Vqt) (q - 1) £ q-( r+2 ^ 2 e-^I r+2j (2^qt), 

3=1 

where I denotes the I-Bessel function. 

Proof. It is immediate that when t = 0, the above series is equal to Sq(t) as required, since 
-£r(0) = for x 7^ and io(0) = 1. Since K(0, r) = <5o(?"); it remains to show that the above series 
is equal to K(t, 0) for t > 0. 

Denote by /(r) = K(t,r) and /(r) = dK(t,r)/dt. If r = 0, then the differential equation (fla)) for 
the heat kernel takes the form 

(« + l)/(0)-(g + l)/(l) + /(0) = 0, 
and for r > the differential equation becomes 

(q + l)/(r) - g/(r + 1) - /(r - 1) + /(r) = 0. 
Let #(r) be the above series expansion times e^ +1 ^, or, when written out, 

oo 

g(r) = q~ r / 2 I r (2^qt) - (q - 1) £ ^T^ 2 ^ 2 /^-^) 

3=1 

It is an elementary exercise to show that the series expansion satisfies the characterizing differential 
equation for the heat kernel if and only if we have the differential equations 

+ 1)5(1) + ff(0) = 0, (6) 

and for every r > 

-qg(r + l)-g(r-l)+g(r)=0. (7) 
Let us first verify ©. We begin by writing the 

oo 

Left-hand-side of © = -(q + l)q- 1 l 2 I 1 {2^/qt) + (q + l)(q - 1) ^ q" (1+2j)/2 /i+ 2j i(2V?*) 

3 = 1 

OO 

+ 2jql' (2jqt) -(q- l)2^qJ2 Q~ 2j ^i 2 ^)- 

3 = 1 
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Using the basic relation I r _i(z) + I r +\(z) = 2I ! r {z) (H), we obtain the expression 

OO 

Left-hand-side of @ = -{q + l)q- 1/2 I 1 (2 y ^t) + (q + l)(q - 1) ^ g~ (1+2j)/2 /i+ 2j <(2jqt) 

+ vWi(V9*) + '-i(V9*)) 

OO 

- (q- 1)^^^(/ 2J+1 (2^) + J 2J ._ 1 (2 V /^)). (8) 

i=i 

Collecting terms, and recalling that i_i = ii, we can evaluate the coefficient of each /-Bessel 
function in (El): 



Jo : there are no Jo terms, 

h ■ -(a + i)<r 1/2 + 2^-(q- i)^- 1 = 0, 

I 2j+1 : (q 2 - l)q-^' 2 - (q - l)Jqq-* - (q - l)^qq~^ 

= (q-l)q- 3 ((q + l)q- 1/2 -q 1/2 -q- 1/2 ) 
= 0. 

In other words, the left-hand-side of © is zero, as required. 

Let us now check the case when r > 0. Again, we begin by writing the 

OO 

Left-hand-side of © = -q l ~ {r+1],2 I r+ i{2^t) + q(q - 1) q' {r+1+2j)/2 I r +i+2j(2^qt) (9) 



-(r 



- 1 )/ 2 / r _ 1 (2Vgt) + (g-i)£?- (r - 1+2j ' )/2 / r -i +2j -(2^) 



J'=l 



2^qq- r /X(2^qt) 2^q(q - 1) £ g -(^')/2j; +2 . (2 ^) 



r'- r / 2 I r+1 (2^t) + q(q - l)^ (? -(''+ 1 + 2 ^/ 2 / r+1+2j (2^) 



- 1 )/ 2 / r _ 1 (2^t) + l)5^g-('- 1 + : W)/ 2 / 1 ._ 1+!b -(2 > /«i) 



J'=l 

OO 

,1/2-r 
-(r 

J'=l 

+ V9<r r/2 a.+i(2^) + 4-i(2V9t)) 

OO 

- y/q(q- l)E^ (r+2j)/2 (4 + l +2j (2^t) +/ r _i +2j (2V^)). (10) 
3=1 

As above, we can evaluate the coefficient of each 7-Bessel function in (J9j> : 

I r - 1 :~q~^/ 2 + ^q- r/2 = 
I r : there are no I r terms, 

I r+1 : -g 1 /2-,/2 + {q l)q -(r + l)/2 + ^ q -r/2 ^ 1)g -(r + 2)/ 2 = Q 

I r+2j+1 : g(g - + {q 1)g -(, + i + 2,)/ 2 

- ^(<Z - l)q-^' 2 Jq{q l) q -lr+V+W 
= {q- l)q-^ 2 q^(q 1 / 2 + q- 1 / 2 - q 1 ' 2 - q- 1 ' 2 ) 
= 0. 

In other words, the left-hand-side of (JJJ is zero, as required, which completes the proof of the 
proposition. □ 

In the following subsection we indicate another approach to the proof of Proposition 13.11 
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3.2 The horospherical transform 

Every geodesic ray 7 in the tree emanating from a fixed base point xq can be viewed as an "ideal 
boundary point at infinity" . To each such 7 there are associated horospheres, one for each integer 
n: 

H„ = {x e X : lim [d(-y(k),x) - k] = n} 

k— >oo 

where d is the natural combinatorial distance in the graph. 

We fix a geodesic ray 7 and may then consider the associated horospherical transform of functions 
/ : X K denoted by 

H/ : Z M 

and defined by H/(n) = X^eH For a ra dial function, decaying fast enough, we have the 

inversion formula 

00 

f(r) = r r (H/)(r) - (« - 1) X; eT (r+2j) (H/)(r + 2j) (11) 

i=i 

for r > 0. This is stated in HNT06 on pages 7-8 for / of finite support. 

If we apply the horospherical transform to the equations (jlal llb[) characterizing the heat kernel, 
we get 

(q + l)HK(t, n) - (qHK(t, n + 1) + HK(t, n-l)) + ^-HK(t, n) = 0, 

at 

for n S Z and with Hif(0, n) — So(n). The solution to this difference-differential equation can be 
seen to be (cf. section 3.2 in |HNT06j or }KN06j ) 

/(i,n)=<r" /2 e- (9+1) %(2^)- 

As already remarked, the heat kernel on a regular tree is radial, so by inserting the above expression 
into the inversion formula (|llj) (here we need to go beyond finitely supported functions for which 
this formula was stated in |HNT06] ) we can get a different proof of Proposition 13. II 

3.3 Heat kernels on regular graphs 

Let q > be an integer and X a (q + l)-regular graph. We fix a base point xq £ VX which we 
will suppress in the notation. The heat kernel on X can be obtained from periodizing the heat 
kernel K q+ i on the universal covering space, the (q + l)-regular tree T q+ i, over the covering group 
r. Following the remarks in subsection 12.31 we have (with a slight abuse of notation) that 

K x (t,x) = Y,K q+ i{t, 1X ). 

Recall that c n (x) denotes the number of paths without backtracking from the identity xq to x of 
length n in X. We also use c° = c n (xo) as notation for the number of geodesic loops, i.e. closed 
paths without backtracking, starting at xq. Note that c n (x) is equal to the number of elements of 
the form 7X for some 7 ST on the radius n sphere in T q+ \. We therefore have 

K x {t,x) = c n (x)K q+ i(t,n) 

n>0 

or more explicitly by inserting the expression from Proposition 13.11 for K q+ i(t,n), 

00 

K x (t, X) = e -^ +1 > X CnW E d q { 3 )q- n/2 -n n+2] {2J-qt), 
n>0 j=0 

where d q (j) is 1 if j = and 1 — q otherwise. A rearrangement of the terms gives 

K x (t,x) = e -^ +1 > 4 X b m (x)q- m / 2 I m (2^qt), 

m>0 
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where b m {x) = c m {x) — (q — l)(c m -2{x) + c m -A{x) + •••) where the last term is c\{x) if m is odd 
and co(x) if m is even. This is also with the understanding that bg(x) = cq(x) and b\{x) = c\(x). 
This proves Thcorem ll.il 

Now specialize to x = xq. In view of Propositions 12. II and 13.11 we obtain Corollary [L2] as well. 
3.4 Spectral theory 

An excellent reference here is that of Mohar and Wocss }MW89| . One has that there are spectral 
measures fi x such that (suppressing xq) 

K x (t,x) = J e- xt dfi x (X). (12) 

In particular if X is a finite graph with n vertices, then the Laplacian has eigenvalues = Ao < 
Ai < ... < A„_i and corresponding orthonormal eigenfunctions <pj. The heat kernel may thus be 
written as 

n— 1 

K x (t, x , x) = -J2 e _Ajt 0j(^)^(^o)- (13) 

71 3=0 

4 Ihara formulas 

In this final section we describe how the G-transform introduced in Section [2.61 applied to the heat 
kernel gives rise to the Ihara zeta function. 



4.1 Zeta functions 

Motivated by Sclbcrg's work, Ihara defined a zeta function for a finite graph X, which is now 
referred to as the Ihara zeta function. The product formula for the Ihara zeta function is 

c^)=n( i - uKP) )" 1 
[p] 

where the product is over equivalence classes of prime geodesies and I the length. (Actually, 
Ihara worked in a specific group setting, but Serre remarked in the preface of |Ser03| that the 
definition could be given a simple interpretation in terms of graphs.) By a general calculation, see 
for example [Terlll p. 29], one has 

00 N 

iogcJ?(«) = E— «™ 

* — ' m 

m—l 

where N rn is the number of closed geodesies of length m. Thus the Ihara function is a zeta type 
function similar to those appearing in the classical works of Artin, Hasse, and Weil on counting 
points of varieties in finite fields. 

The numbers N m are not defined for infinite graphs X. In the case of transitive graphs, a natural 
replacement is iV^, since in the finite case one has N m = nN® n , where n is the number of vertices. 
So we can define a zeta function for any (not necessarily finite) vertex transitive graph via 

m— 1 

More exotically one could define a two variable function f via 

b m (x) 



log Of (u, a;) = ^2 



-u 



Following Riemann one has another set of zeta functions by instead taking the Mellin transform 
of the the heat kernel. This is a subject for another paper. 
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4.2 The periodization side 

We will apply the transform 

/>oo 

Gf(u) = (u- 2 -q) er^+^^e^+^fitjdt 
Jo 

first to our heat kernel expression 

oo 

K x (t,x ,x) = e-^+D* b m (x)q- m / 2 I m (2^t) 

of Theorem 11.11 Using the basic formula in subsection I2.6[ in the case x xo, so bo(x) = we 
have that the G-transform of the heat kernel is equal to 

- £ b m (x)u m = A f] = -?-lo g ( x (u,x). 

u ^— ' ou m ou 

rn—O m— 1 

In the case x = xq, the G-transform of the heat kernel on the diagonal is equal to 



oo 



u z — ' ou z — ' m ou 

m— m— 1 

In the vertex transitive case, for the case x = Xq we apply the transform to the expression of 
Corollary O 

OO 

K x (t,x ,x ) = K q+1 (t,x ) + e-(« +1 >* K n q~ m / 2 I m (2^t) 

771—1 

which gives (summing the geometric series arising from the first term) 

1 oo 
rn—1 

This can in turn be written as 

#- log« + ^f lo g (l-^) + A£^ 
au 2 ou ou * — ' m 

m— 1 

= ^ io gu + — ^ mi - . 2 ) + - io g (aw). 

We have proven 

Proposition 4.1. For X a (q + l)-regular vertex transitive graph, 



(GK X )(;X )( U ) = -^- 

ou 



logu+^— log(l - u 2 ) + \og( x (u). 



In summary, the G-transform of the heat kernel yields expressions involving the Ihara zcta function 
together with other trivial terms. In the setting of compact quotients of rank one symmetric spaces, 
there is a similar change of variables in the Laplace transform so that when applied to the trace 
of the heat kernel, one obtains the Selberg zeta function. With this in mind, our approach to the 
Ihara zeta function as an integral transform of the heat kernel is in line with a known approach 
to the Selberg zeta function in many settings. 
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4.3 Ihara's determinantal formula 

Now wc deduce the classical Ihara determinantal formula. Let X be a (g + l)-rcgular graph of 
finite vertex cardinality n. From (|13p . 



iCx(t,a;o) = - V"e 

r7 * — ' 



n 



The G-transform of the righthand side is a simple integration which yields 
1 n ~ 1 l 1 r) n— 1 1 

V-9)E— Tr4TrTT = - i ^E 1 °g-( 1 -(« + 1 - A >+^)- 

n qu + 1/u — (q + 1 — A, ) n ou m v ' 

Comparing t his last expression with Proposition ^. 1 1 from the pcriodization side (and also verifying 
that integration constants match up) we immediately get Ihara's formula, namely 



since, as remarked above, 



= (1 - u 2 )«(9-i)/2 dct((l - (q + l)u + qu 2 )I + Au), 



m— 1 



This formula is also known to hold more generally for non-regular graphs, see the references 
mentioned in the introduction. 



4.4 First extension of Ihara's formula 



Using spectral theory, we obtain a similar formula for infinite transitive graphs, this time using 
our zeta function instead of Ihara's. The common point is that both zeta functions are obtained 
as G-transforms of the heat kernel, and the determinantal formula follows from having another 
expression for the heat kernel, namely that which comes from spectral theory. 
We use the notation fi = fi XQ . Since the functions involved are positive, we may change the order 
of integration in our integral transforms and arrive at the expression Equating the expression in 
Proposition 14.11 with the G-transform of the spectral expansion of the heat kernel given in (fT^j) . 
we arrive at 



d_ 

du 



log it ■ 



(I 



1 



log(l - it 2 ) +logCx(u) 



— / log - (1 - (q + 1 - X)u + qu 2 ) dfx(\). 
ou J u 



Note: since the functions involved are positive, we are justified in interchanging the spectral 
integral with the G-transform integral on the righthand side. We now integrate this equality, 
noting that at u — both sides are to determine the integration constants. We get the formula 

logw + log(l - u 2 ) + logCx(w) = - /log i (1 - (g+ 1 - \)u+ qu 2 ) rf/i(A), 



which leads to 



cxiuy^ii-u 2 )^ 2 



exp 



J log (l - (q + 1 - X)u + qu 2 ) dn(X) 



(14) 



This is Theorem 1 1.31 We further remark that equation ([TJJ clearly generalizes the Ihara determi- 
nant formula since for vertex transitive graphs with a finite number n vertices one has fj^ = (3c • 
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4.5 Second extension of Ihara's formula 



Here we do not specialize to x = xq. The resulting identities involve counting geodesies paths, 
not only closed geodesies paths. Alternatively, as in the most classical situation, our considera- 
tion corresponds to computing the Hurwitz zeta function instead of the Riemann zeta function. 
With the same calculations as above, one gets in the finite graph case, (at one point one uses 
orthogonality of eigenfunctions) the formula 



Thus the eigenfunctions come in to determine the more precise count of geodesies. The lead 
asymptotic as the length goes to infinity behaves the same as for the closed geodesies since the 
trivial eigenvalue has the constant function as eigenfunction. From an intuitive viewpoint, this 
observation is clear: For a fixed x and large length m the geodesies do not look much different 
from a closed geodesic. In symbols, if m 3> 1 then i«io. 

We have the analogous formula for infinite regular graphs, namely that 



References 

[Ahu87] Guido Ahumada. Fonctions pcriodiqucs ct formule des traces de Sclberg sur les arbres. 
C. R. Acad. Sci. Paris Ser. I Math., 305(16):709-712, 1987. 

[Bas92] Hyman Bass. The Ihara-Selberg zeta function of a tree lattice. Internat. J. Math., 
3(6):717-797, 1992. 

[Bcd97] Debe Bcdnarchak. Heat kernel for regular trees. In Harmonic functions on trees and 
buildings (New York, 1995), volume 206 of Contemp. Math., pages 111-112. Amer. Math. 
Soc, Providence, RI, 1997. 

[CJK10] Gautam Chinta, Jay Jorgenson, and Anders Karlsson. Zeta functions, heat kernels, 
and spectral asymptotics on degenerating families of discrete tori. Nagoya Math. J., 
198:121-172, 2010. 

[CMS00] Michael Cowling, Stefano Meda, and Alberto G. Setti. Estimates for functions of the 
Laplace operator on homogeneous trees. Trans. Amer. Math. Soc, 352(9):4271-4293, 



[CMS01] Bryan Clair and Shahriar Mokhtari-Sharghi. Zeta functions of discrete groups acting on 
trees. J. Algebra, 237(2):591-620, 2001. 

[CY99] Fan Chung and S.-T. Yau. Coverings, heat kernels and spanning trees. Electron. J. 
Combin., 6:Rescarch Paper 12, 21 pp. (electronic), 1999. 

[DM06] Jozcf Dodziuk and Varghese Mathai. Kato's inequality and asymptotic spectral prop- 
erties for discrete magnetic Laplacians. In The ubiquitous heat kernel, volume 398 of 
Contemp. Math., pages 69-81. Amer. Math. Soc, Providence, RI, 2006. 

[Dod06] Jozcf Dodziuk. Elliptic operators on infinite graphs. In Analysis, geometry and topology 
of elliptic operators, pages 353-368. World Sci. Publ., Hackcnsack, NJ, 2006. 

[FZ99] Dominique Foata and Doron Zeilberger. A combinatorial proof of Bass's evaluations of 
the Ihara-Selberg zeta function for graphs. Trans. Amer. Math. Soc, 351(6):2257-2274, 



\og(x(u,x) = ~Y, /jW/j^o) iogl 1 - (g + 1 - A,-)« + qu 2 ). 




2000. 



1999. 



13 



[GIL08] Daniele Guido, Tommaso Isola, and Michel L. Lapidus. Ihara's zeta function for periodic 
graphs and its approximation in the amenable case. J. Fund. Anal., 255(6):1339-1361, 
2008. 

[GZ04] Rostislav I. Grigorchuk and Andrzej Zuk. The Ihara zeta function of infinite graphs, 
the KNS spectral measure and integrable maps. In Random walks and geometry, pages 
141-180. Walter de Gruyter GmbH & Co. KG, Berlin, 2004. 

[Has89] Ki-ichiro Hashimoto. Zeta functions of finite graphs and representations of p-adic groups. 

In Automorphic forms and geometry of arithmetic varieties, volume 15 of Adv. Stud. Pure 
Math., pages 211-280. Academic Press, Boston, MA, 1989. 

[HNT06] Matthew D. Horton, Derek B. Newland, and Audrey A. Terras. The contest between the 
kernels in the Selberg trace formula for the {q + l)-regular tree. In The ubiquitous heat 
kernel, volume 398 of Contemp. Math., pages 265-293. Amcr. Math. Soc, Providence, 
RI, 2006. 

[Iha66] Yasutaka Ihara. On discrete subgroups of the two by two projective linear group over 
p-adic fields. J. Math. Soc. Japan, 18:219-235, 1966. 

[JL01] Jay Jorgcnson and Serge Lang. The ubiquitous heat kernel. In Mathematics unlimited — 
2001 and beyond, pages 655-683. Springer, Berlin, 2001. 

[KN06] Anders Karlsson and Markus Neuhauser. Heat kernels, theta identities, and zeta func- 
tions on cyclic groups. In Topological and asymptotic aspects of group theory, volume 
394 of Contemp. Math., pages 177-189. Amer. Math. Soc, Providence, RI, 2006. 

[KS00] Motoko Kotani and Toshikazu Sunada. Zeta functions of finite graphs. J. Math. Sci. 
Univ. Tokyo, 7(l):7-25, 2000. 

[McK72] H. P. McKean. Selberg's trace formula as applied to a compact Riemann surface. Comm. 
Pure Appl. Math., 25:225-246, 1972. 

[Mne07] P. Mnev. Discrete path integral approach to the Selberg trace formula for regular graphs. 
Comm. Math. Phys., 274(1):233-241, 2007. 

[MW89] Bojan Mohar and Wolfgang Wocss. A survey on spectra of infinite graphs. Bull. London 
Math. Soc, 21(3):209-234, 1989. 

[OB73] Fritz Oberhettinger and Larry Badii. Tables of Laplace transforms. Springer- Verlag, 
New York, 1973. 

[Pal99] B. V. Pal'tscv. On two-sided estimates, uniform with respect to the real argument and 
index, for modified Bcsscl functions. Mat. Zametki, 65(5):681-692, 1999. 

[Ser97] Jean-Pierre Serre. Repartition asymptotique des valeurs propres de l'operateur de Hecke 
T p . J. Amer. Math. Soc, 10(1):75-102, 1997. 

[Ser03] Jean-Pierre Serre. Trees. Springer Monographs in Mathematics. Springer- Verlag, Berlin, 
2003. Translated from the French original by John Stillwell, Corrected 2nd printing of 
the 1980 English translation. 

[Sun86] Toshikazu Sunada. L-functions in geometry and some applications. In Curvature and 
topology of Riemannian manifolds (Katata, 1985), volume 1201 of Lecture Notes in 
Math., pages 266-284. Springer, Berlin, 1986. 

[Sun94] Toshikazu Sunada. Fundamental groups and Laplacians [ MR0922018 (89d:58128)]. In 
Selected papers on number theory, algebraic geometry, and differential geometry, volume 
160 of Amer. Math. Soc. Transl. Ser. 2, pages 19-32. Amer. Math. Soc, Providence, RI, 
1994. 



14 



[Sun08] Toshikazu Sunada. Discrete geometric analysis. In Analysis on graphs and its ap- 
plications, volume 77 of Proc. Sympos. Pure Math., pages 51-83. Amer. Math. Soc, 
Providence, RJ, 2008. 

[Tcrll] Audrey Terras. Z eta functions of graphs, volume 128 of Cambridge Studies in Advanced 
Mathematics. Cambridge University Press, Cambridge, 2011. A stroll through the 
garden. 

[TS07] A. A. Terras and H. M. Stark. Zcta functions of finite graphs and coverings. III. Adv. 
Math., 208(l):467-489, 2007. 

[TW03] Audrey Terras and Dorothy Wallace. Selberg's trace formula on the fc-regular tree and 
applications. Int. J. Math. Math. ScL, (8):501-526, 2003. 



Gautam C hint a 
Department of Mathematics 
The City College of New York 
Convent Avenue at 138th Street 
New York, NY 10031 U.S.A. 
e-mail: chinta@sci.ccny.cuny.edu 

Jay Jorgenson 
Department of Mathematics 
The City College of New York 
Convent Avenue at 138th Street 
New York, NY 10031 U.S.A. 
e-mail: jjorgenson@mindspring.com 



Anders Karlsson 
Mathematics Department 
University of Geneva 
1211 Geneva, Switzerland 
e-mail: anders.karlsson@unige.ch 



15 



